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Abstract 

We find an asymptotic expression for the first eigenvalue of the 
biharmonic operator on a long thin rectangle. This is done by find- 
ing lower and upper bounds which become increasingly accurate with 
increasing length. The lower bound is found by algebraic manipula- 
tion of the operator, and the upper bound is found by minimising the 
quadratic form for the operator over a test space consisting of sepa- 
rable functions. These bounds can be used to show that the negative 
part of the groundstate is small. 

1. Introduction 

There is a considerable literature on studies of the biharmonic operator act- 
ing in L 2 (Q) for particular regions Q C IR 2 , especially the square, disk and 
punctured disk. In this paper we study the operator acting in L 2 (Rh), where 
Rh is the rectangle [0, h] x [0,1]. Difficulties in studying the biharmonic 
operator arise because the eigenvalue problem 

A 2 / = fif, f = |J = on dR h (1) 

is not exactly soluble. The boundary conditions in this problem are called 
Dirichlet boundary conditions or clamped plate boundary conditions. Nu- 
merical analysts have succeeded in proving a number of interesting results 
about the groundstate of the biharmonic operator and the corresponding 
eigenvalue for the square and some rectangles, however there are very few 



1 



previous results concerning the /i-dependence of spectral quantities. One 
such result is by Behnke and Mertins [|BM95|| . See note [|. 

For the unit square the best current enclosure 



\i x = 1294.9339 
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for the first eigenvalue is due to C. Wieners [ Wie using the Tempel-Lehmann- 
Gorisch method to obtain the lower bound, and minimisation of the quadratic 
form of the operator on a certain space of test functions for the upper bound. 
The enclosure is guaranteed by interval arithmetic programming. The accu- 
racy with which this value has been computed has increased with computing 
power over the last sixty years. In 1937 Weinstein ||Wei37|| introduced a 



method which theoretically enabled him to calculate a limiting sequence of 
lower bounds for yU 1; although his hand calculations 

1294.956 < /ii < 1302.360 

were slightly inaccurate. Results of intermediate accuracy have been given 
by Aronszajn |[Aro51|| (from a citation in [Fic78| , page 74]), Bazley, Fox and 
Stadter fBFS67j, De Vito, Fichera, Fusciardi and Scharf |VFFS66|1 and many 
others. 

The biharmonic operator, which we shall denote by A 2 |dir, is defined as the 
non-negative self-adjoint operator associated with the closed quadratic form 



Q(f) 



7 |A/| 2 , if/eiy 2 ' 2 (i4); 



(2) 



oo, 



otherwise. 



See |pav95b| , theorem 4.4.2] for details. Using the Rayleigh-Ritz formula ||Dav95b 
section 4.5], the first eigenvalue of the biharmonic operator is given by the 
expression 



^{h) = inf{Q(/) : f E L 2 (R h 



!}• 



(3) 



We give formulae for lower and upper bounds \i(h), Vi(h) for fii(h) and use 
their asymptotic expressions to prove that 



^(h) = c 4 + 2dc 2 7r 2 h- 2 + 0(h- s ) 
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as h — > oo, where c ~ 4.73004 is the first positive solution of the transcen- 
dental equation 

cosh c cose = 1, 

and (4) 

2 tanh c tan c — c tanh c — c tan c 

d = w 0.54988. 

ctanhc — ctanc 

Elementary algebraic manipulation is used to find the lower bound; no benefit 
is derived for large h by using more involved methods such as Weinstein's 
truncation of operators. The upper bound is found, as in most other papers, 
by minimisation of the quadratic form of the operator over a certain test space 
of functions. Observing that eigenfunctions of the biharmonic operator are 
close to being separable functions, we choose our test space to consist of all 
such functions. This simplicity of approach allows us to find the asymptotic 
formulae. The lower and upper bounds we find are also useful for small values 
of h. They are within 0.72% of the actual eigenvalue for all h G [1, oo). See 
figure |3| and table |l|. 

The advantage of our results lies in the fact that they are valid for all values of 
h, and also in the limit as h —>■ oo. They are also simple to compute, without 
the need for finite element algorithms. The natural traded disadvantage 
is that for a particular rectangle with comparable side lengths, where h is 
close to 1, numerical analysts are able to use powerful computers using finite 
element methods to compute eigenvalues to a far higher degree of accuracy. 

A fundamental issue in the study of fourth order operators is the fact that 
the corresponding semigroup is not positivity preserving. This is exhibited 
by non-positivity of the groundstate of the biharmonic operator for certain 
domains Q, a feature first noticed by Bauer and Reiss (1972) [ |BR72j| for the 
square. A rigorous proof of this fact has been given by Wieners (1995) |[Wic96 



who finds a function which is slightly negative very near the corners of the 
square and pointwise extremely close to the actual groundstate. Kozlov, 
Kondrat'ev, and Maz'ya (1990) | |KK1V190[ had in fact already obtained a 



more informative result. They managed to show that if the region has an 
internal angle of less than 146.30° then the groundstate oscillates infinitely 
often in sign as one approaches the corner, although there are no explicit 
data concerning where the first oscillation occurs. 
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Using the bounds A 1; V\ and also a lower bound A 3 on the third eigenvalue, 
we prove the bound 

jjfn.<£z^ = 0(fc -v,) (5) 

||/l 1 1 25 (A 3 - V\) VI1 

as /i — > 00, on the size of the negative part /f of the groundstate /1 of A 2 |dir. 
See figure £| and table for a plot and values of this function for small values 
of h. In particular we see that for the case of the square, 

l^J- 2 < 0.0484. 
I l/i II2 

By 'negative part of the groundstate' we mean the negative part when the 
eigenfunction is positive in the centre of the square. 

It is possible to use the L 2 bound and a Sobolev embedding theorem to imply 
that 

II fi II 00 < (?A - Ai) 1/2 A3 4 _ n(h -i/2\ 
Wkh - 2(A 3 -^ - U[h } 

as h — > 00, and so the negative part of the eigenfunction /1, already small for 
h — 1, vanishes asymptotically as h — > 00 in the L 2 and L°° senses. However 
since 

II f ,11. _ 

0(/T 



II A II 



IIAII2 

only the L 2 result (S) is of particular value. See note [L5| for details. 



2. A Fourth Order Operator in one Dimension 

Our approach to analysis of the biharmonic operator acting in L 2 (Rh) in- 
volves attempting to separate variables. Eigenf unctions of A 2 |dir are not 
separable functions, a fact confirmed by oscillations at the corners, so it is 
remarkable that we obtain such good estimates for the first eigenvalue. The 
method is successful because the eigenfunctions are close to being separa- 
ble. The following sections will rely heavily upon spectral analysis of the 
self-adjoint operator 
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acting in L 2 ([0, h]), with quadratic form domain Wo' 2 ([0, h]). Since H(h,a) 
is bounded below and has compact resolvent we may order the eigenvalues 
as an increasing list. 

Let a(h, a, n) : (0, oo) x R x N — > K. be the function which associates to each 
a the n-th eigenvalue of the operator H(h,a) acting in L 2 ([0,h}), and let 
p n (a) = a(l,a,n). 

Theorem 1: 

(i) The first eigenvalue p\ of H(l, a) is an increasing and concave function 
of a. The functions a and p are related by the equation 



For a > 0, pi is analytic; 

(ii) For a > the Green's function of H(h, a) is positive. It follows that 
for such a the first eigenvalue is of multiplicity one and the groundstate 
is positive; 

(hi) For a > let f be the n-th eigenfunction of H(h, a) and have unit L 2 
norm. Then 



(iv) For a > let (3 > a be the n-th solution of the transcendental equation 




(7) 



\\f'\\l = \p'n(h 2 a)h- 2 ; 



(8) 



cosh y/ p + a cos \ ffi 



a 



sinh y/ P + a sin 



— a — 



— a = 1. 



— a 



(9) 



Then p n (a) = j3 2 — a 



2. 



(v) The following are asymptotic expansions of p n at and oo: 



p n {a) = 2n 2 n 2 a + 4 v / 2n 2 7r 2 a 1/2 + 0(1) 
p' n (a) = 2n 2 vr 2 + 2 v / 2nVoT 1/2 + C^oT 1 ) 



as a 



as a 



oo: 



oo: 




as a -> 0; (10) 



as a 



0; 
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Figure 1: Plot of the first four eigenvalues of H(l, a) 

where c n is the n-th positive solution of the equation cosh c cose = 1, 
and 

2 tanh c n tan c n — c n tanh c n — c n tan c n 
c n tanh c n c n tan c n 

(vi) The function p[ is monotone decreasing and 

2tt 2 < p\(a) < 2d lC \ (11) 

for all a > 0. 

Figure [I] shows the first four eigenvalues of H(l, a). It was plotted by Mathe- 
matics, using the implicit formula (||) in theorem [I] (iv) for positive p, with 
a similar formula for negative p. 

Note 2: The portion of the graph for negative a is not relevant in this 
paper but has been included to demonstrate that even simple fourth order 
operators have completely different eigenvalue behaviour to the second order 
theory. Coincidence of p-m-i an d P2n occurs for 

(a, p) = (-(m 2 + n 2 )7r 2 , -(m 2 - n 2 ) 2 vr 4 ), (12) 

where m — n G N. □ 

Proof of theorem [l] (i): The function p\ is increasing because the perturb- 
ing operator is positive. Using the Rayleigh-Ritz formula, 

Pi(a) = mi{{H(l,a)g,g) 1 : \\g\\ 2 = l,ge Dom(if(l, a))}. 
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Suppose that a = A«i + (1 — \)a where < A < 1. Then 

(H(l,a)g,g) 1 = X(H(l,a 1 )g,g) 1 + (1 - X){H{l ) a Q )g ) g) l . 

Let e > and choose g G Dom(if (1, a)) such that \\g\\2 = 1 and (H(l, a)g, g)\ < 
Pi(a) + e. Then 

Pi(a) + e > X(H(l,a 1 )g,g) 1 + (1 - A)(#(l, a )g,g)i 
> Api(ai) + (1 - A)pi(a ). 

Since this holds for all positive e, is concave. 

For fixed a let /„ G L 2 ([0,/i]) be the n-th eigenfunction of H(h,a). Define 
g n G L 2 ([0, 1]) by # n (y) = f n (hy). By the chain rule, 

H(l, h 2 a)g n = - 2h 2 a^^ = h*a(h, a, n)g n , (13) 

so g n is an eigenfunction of H(l,h 2 a) with eigenvalue h A a(h,a,n). It fol- 
lows that a 4 cr(a, a, n) > p n (h 2 a). By a similar reverse argument we obtain 
equality. 

The family H(l, a) of differential operators indexed by a is a holomorphic 
family of type B and so p n are analytic except where the eigenvalues swap. 
See [|Kat80| , chapter VII §4] . We shall see in part (ii) that swapping does not 
occur when a is positive. □ 

Proof of theorem [l] (ii): For a > define G : [0, l] 2 — > K by 

G(x,y) = <^ (14) 
[-A;(x, 1 -y), x<y 

where 

c = a 3 (2(l — cosh a) + a sinh a) (15) 

and 

k(x, y) = (sinh a — a) (cosh ax — 1) (cosh ay — 1) 

— (cosh a — 1) (cosh ax — 1) (sinh ay — ay) 

— (cosh a — 1) (sinh ax — ax) (cosh ay — 1) 
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+ sinha(sinhax — ax)(sinh ay — ay). (16) 
For g E L 2 ([0,l]) define / by 



Putting 



the identities 



f( x ) = / G(x,y)g(y)dy. 
Jo 



k(l-x,x)- k(x,l-x) = k^°\x,y) = a 2 k^°\x,y) 
-k^°\l-x,x)-k^°\x,l-x) = 

k^°\l-x,x)-k {2 '°\x,l-x) = k{0,y) = 

-k^°\l-x,x)-k {3 '°\x,l-x) = c k^' o \0,y) = (17) 

imply that 

/(0) = f'(0) = /(l) = f (1) = 

and 

dx 4 dx 2 ' 
so G is the Green's function of the operator 

d 4 2 d 2 
- a 2 



dx A dx 2 

acting in L 2 ([0, 1]) with Dirichlet boundary conditions. For positivity of 
G(x,y) we need to to prove that k(x,y) is positive in the triangular region 

T = {(x,y) : x,y > and x + y < 1}, (18) 

In order to do this, we introduce the function 

, / x , 1 ( sinha; — x\ 



The reader should verify that 



lim d)'(x) = - 



and 

fx/2 



r-x/2 ■ I ( r \ 

x) — I 4 sinh 2y [sinh 2 ytanhy — y 3 ]dy / I / y(coshy — l)dy I >0 



for x > 0. These features imply that is convex and increasing for positive 
x. Let (x, y) G T. Then 

(f)(a) > (f)(ax) + <f)(ay). (20) 
Hence using inequality (EDf) and a two angle tanh identity, 



Thus 



sinh a — a sinh a — a sinh ax — ax sinh ay — ay 
cosh a — 1 cosh a — 1 cosh ax — 1 cosh ay — 1 
= tanh{0(a)}[l + tanh{0(aa;)} tanh{0(a?/)}] 
> tanh{0(ax) + (p(ay)}[l + tanh{0(ax)} tanh{0(ay)}] 
= tanh{0(ax)} + tanh{0(ay)} 
sinh ax — ax sinh ay — ay 
cosh ax — 1 cosh ay — 1 

k(x,y) = [(sinh a — a) (cosh ax — 1) (cosh ay — 1) 

— (cosh a — 1) (cosh ax — 1) (sinh ay — ay) 

— (cosh a — 1 ) (sinh ax — ax) (cosh ay — 1 ) 
+ (sinh a — a) (sinh ax — ax) (sinh ay — ay)] 
+ a(sinh ax — ax) (sinh ay — ay) > 

as required. 

It follows that the Green's function of H(h, a) for all h, a > is positive 
because of the relationship (|i~3|) between if(a, a) and H(l, h 2 a) established 
in part (i). 

See [|RS78| , theorem XIII. 44] for a proof that the groundstate of H (h, a) is 
positive and the associated eigenvalue is multiplicity one. □ 

Proof of theorem I (iii): For a > let / = f a e Wq' 2 ([0, h])nC°° satisfy 

d 4 / d 2 f 

dx 1 ~ 2a dx^ = a ( h,a,n ' )f (21) 
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and have unit L 2 -norm. Then f a is a critical point of the functional 



>" 2a <d^ 



dx 2 ' dx 2 



in the sense that if ip(t) is C l with respect to £, 
i[>(0) = f a then 



(22) 

1 for all t and 



0. 



t=o 



Letting = f a+t and using the relationship (|7|) between p and a, 

p' n {h 2 a)h~ 2 = —a(h, a,n) 
da 

= -T-£a(fa) 



da 
d 



dt 



t=a 



t=0 



□ 



Proof of theorem |T| (iv): Since H(l,a) is positive for a positive, eigen- 
functions may be found by solving the auxiliary equation y A — 2ay 2 — p = 
where p > 0. Since a — \J a 2 + p < < a + \/ a 2 + p, there are four distinct 
roots, two real and two imaginary. Denoting these roots by a, —a, ib, —ib, we 
see that a 2 — b 2 = 2a and a 2 b 2 = p. 

There exists a combination of the functions cosh ax, sinh ax, cos bx, sin bx 
which satisfy the boundary conditions of H(l, a) if and only if 



det 



/ 1 1 \ 

a b 

cosh a sinh a cos b sin b 

ya sinh a a cosh a —b sin b b cos b ) 



0. 



Simplifying this determinant, the equation becomes 
2ab cosh a cos b — (a 2 — 6 2 ) sinh a sin b 



2ab, 
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so p = f3 2 — a 2 is an eigenvalue of the operator if and only if 



cosh \f]3 + a cos y/ /3 — a = sinh \J (3 + a sin y/ j3 — a = 1. 

\J (3 2 — a 2 



□ 



Proof of theorem [1] (v): The proof of asymptotic formulae for p and p' 
as a — > oo is given in the appendix. It is possible to find the asymptotic 
formula for p as a — > by a similar method. The formula for p' then follows 
by differentiation because p is analytic at 0. Here we give a sketch of an 
alternative proof of the case a — > for interest: 

Let 7„ = cosh c n — cos c n and <5 n = sinh c n — sin c n . Define 

g n ( x ) — c °sh c„x — sinh c n x — cos c n x + ^ sin c n x. (23) 

We claim that (<7 n )neN is an orthonormal sequence of eigenfunctions of H(l, 0), 
and so p n (0) = c*. These eigenfunctions represent the fundamental modes 
of the clamped beam. We claim that, \\g' n \\ 2 = dnC^, and so from (§) we see 
that p' n {0) = 2d n c 2 l . 

Verification of each of these claims is not trivial; indeed a lengthy calculation 
is needed even to establish that \\g n \\2 — 1 for each n. This task is left to the 
reader. □ 

Proof of theorem |l] (vi): Since p\ is concave, p' x is decreasing. The result 
now follows immediately from part (v). □ 



3. Lower Bounds on Eigenvalues 

We find lower bounds \i(h), A 3 (/i) for pi(h) , p 3 (h) respectively, by elementary 
algebraic manipulation of the biharmonic operator. This method would be 
referred to as a finite renormalisation procedure in the physics literature. 

Theorem 3: 

Ai(/i):= Pl (7r 2 h 2 )h- 4 + Pl (7r 2 h- 2 ) -2vr 4 ^ 2 < p^h), (24) 

WM . ( Pl (n 2 h 2 )h-* + p 2 (n 2 h- 2 )-2n*h- 2 ,\ 
A3 ^ :=mm (p 3 (vr^)^ + p 1 (vr^)-2vr^ \ ~ ^ ^ 

where p n are defined in theorem [J. 
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Proof: In this proof we consider, where relevant, restrictions of operators to 
C^°([0,h}). Let A h denote the biharmonic operator acting in L 2 ([0,h}), and 
let Bh denote the Dirichlet Laplacian acting in L 2 ([0, h]). Then 

A 2 = A h <g> lj + l h <g> A x + 25^ ® fix 

= + 2-K 2 B h ) ® li + l fc ® (Ai + 2vr 2 /i- 2 5 1 ) 

+ 2(5 ft - 7r 2 /r 2 l ft ) ® (Si - tt 2 1i) - 27r 4 /T 2 l h ® li- 

The operator — 7i 2 h~ 2 l h ) ® (5 X — 7r 2 lx) has eigenvalues 

vr 4 /i- 2 (m 2 - l)(n 2 - 1) m, n G N, 

with the corresponding complete orthonormal sequence of eigenfunctions 

2h~ 1 ^ 2 sin rrmhT x x sin tot?/, 

and so is non-negative. Hence 

A 2 | Dm > {A h + 2ix 2 B h ) ® li + 1 A <g> (Ax + 2vr 2 / i - 2 5 1 ) - 27r 4 /T 2 l h <g> 1l 

It now follows that for large enough 

»n(h) > <r(h, 7T 2 , n) + tr(l, tt 2 /*" 2 , 1) - 2tt 4 ^ 2 
= p n (vr 2 /i 2 )/r 4 + pi(vr 2 /r 2 ) - 2vr 4 /T 2 
= A n (/i). □ 

The functions X n may be explicitly calculated using formulae @, (24) and 
(25). See figures ||, |3] and table [I]. 

Note 4: Despite the fact that the lower bounds A n cross each other, swapping 
of eigenvalues is a not actually a genuine feature of the increasing size of the 
rectangle. In |[BM95|| , Behnke and Mertins show that the eigenvalues veer 



away from each other just before the points where one might expect them to 
cross. □ 



4. Upper Bound on the First Eigenvalue 

We shall find an upper bound v\{h) for fj,\(h) by approximation of the ground- 
state eigenfunction with separable functions. This is called the Hartree 
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Figure 2: Graph of X\(h) and Xs(h) 



method in the physics/chemistry literature and was used in [ Pav95a|| to es- 
timate the groundstate energy of a somewhat similar but simpler problem. 
Define the functional 



£:L 2 ([0,h])xL 2 ([0,l}) 



by 

£{f,g) = Q(f®g) 



h pi 




\A(f(x)g(y))\ 2 dxdy, f e W 2 > 2 ([0, ^])and g e W^([0, 1]); 



r2,2, 



otherwise. 



(26) 



Let 



v^h) = w£{£{f,g) : / E L 2 ([0,h}),g G L 2 ([0, 1]), 



M| 2 = i}, 



(27) 



be our upper bound on fii(h), where the two norms are taken in L 2 ([0,h}) 
and L 2 ([0, 1]) respectively. The remainder of this section is devoted to giving 
a characterisation of u±, and finally an implicit formula. 

Lemma 5: The infimum in the expression ( |27[ ) for V\ is attained. 

Proof: Let Ah denote the biharmonic operator acting in L 2 ([0,h]). Since 
Ah has compact resolvent, it has a complete orthonormal sequence of eigen- 
functions whose corresponding eigenvalues form a divergent sequence. Using 
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this orthonormal sequence we see that the set 

S h :={fEL 2 ([0,h}):\\f\\l + Q h (f)<l} 

is compact. The set 

S :={/ G L 2 ([0, h)) : ||/|| 2 + Q h (f) < 1, > 1/(21* + 1)} 
xfee L 2 ([0, 1]) : \\g\\ 2 2 + Q 1 (g) < 1, \\g\\ 2 2 > 1/(2^ + 1)} 

(28) 

is a closed subset of S(, x 5i C L 2 (R h ), separated from the origin, and so is 
a compact subset of L 2 (R h ) \ {0}. Since the map 

/"W)/ll/lla (29) 

is lower semicontinuous on L 2 (Rh)\{0} it attains its infimum when restricted 
to S. 

Suppose that S(f,g) < 2v\ where / and g have unit norm. We may rescale 
/ and g so that ||/||| + Q h (f) = 1 and \\g\\ 2 + Q 1 (g) = 1. Now 

I ~ \\f\\ 2 2 _ Q h (f) ^ S(f,g) ^ 

11/112 11/112 " II/II2IMII - u 

so H/H2 > l/(2z/i + 1). A similar argument for g shows that (/, g) G S. It 
follows that 

min ^^ii :(/ ' 5)G5 }^ l(/l) ' 



Let (/, g) G S take the minimum value of the map (|29| ) restricted to S. We 
may rescale so that ||/||2 = ||g||2 = 1- Now 

**(/*) < £(f,g) = min ( „ly^ 2 : if, 9) e ^} < n 



From this point onwards, we shall assume that £ is only applied to functions 
of unit norm. We may rewrite the formula for £ as 

{J,9> W'dx 2 ^ + W'^dy 2 '^ 1 W'dy 2 ' 1 ' (3Q) 
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where (., .)h denotes the inner product on L 2 ([0, h]). We shall identify / and 
g by searching for the critical points of S. 

Lemma 6: Let f G Wg' 2 ([0, h\) and g G VFo' 2 ([0, 1]) minimise £(f, g). Then 
f G C°°([0,h]) fl VFo' 2 ([0, h]) is the groundstate of the operator H(h,a g ), 
where 

a g := \\g'\\ 2 2 . (31) 



Proof: Let / G C c °°([0, h\) be such that (/, f) h = 0. Define 

/ + */ 



/(*) == 



By differentiating we see that 



n/+*/ir 



(32) 



/(0) = / and ^(0) = /. 
The minimum of S will be a critical point so 



o=|(/W, 9 ) 



t=0 



9R /d 4 / _ d 2 / . 



'dx 2: 



where the fourth derivative has been taken in the distributional sense. Re- 
placing / by if we see that 



for all such /. It now follows that 

d 4 / . d 2 / 

2a n ^r = A*/ 



(33) 



dx 4 



'dx 2 



(34) 



For suppose otherwise, then there exist /i,/2 G with (f,fi)h and 
(/, /2)h 7^ 0, and /i X ^ /i 2 such that 



dx 2 



1,2. 
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Now let / = (/, f 2 ) h h - (f, fx) h f 2 . Then 

(/, f)h = 

and 

d 4 f d 2 f ~ 

= (^i ~ ^){fJi)h{f,h)h 

7^0, 

which contradicts (|33|). 

It follows from equation ( 151 ) that / G PU m ' 2 for all m G N, and consequently 
/ is smooth. From ( |30D the value of £ at the critical point is 

In order for this to be the minimum, / must be the groundstate of H(h, a g ). □ 

The following theorem gives a formula which characterises the abstract defi- 
nition (|2?D of v\. 

Theorem 7: Let V\ be the upper bound on fix defined by ft2/[) . Then 



//,(//) = f>y(Ir<y)lr l + ( ipi(/i 2 a)/i- 2 ) -// L (/rn)/r-n. (35) 



where a G [7T , dc ] is a solution of the equation 



2a = p' 1 (^p[(h 2 a)h- 2 y 



(36) 



Proof: Let (/, g) be a pair of functions which attain the minimum of £. 
Then by lemma || / is the groundstate of the operator H(h, a g ). Using (g) 
we see that 



By an identical argument, 



«/ = \p',{h 2 a g )h-\ (37) 



a 9 = \p'i{ a j)i (38) 
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1.0 


1286.66 


1295.93 


0.720 


1.2 


940.070 


946.421 


0.676 


1.4 


776.088 


780.618 


0.584 


1.6 


687.796 


691.129 


0.485 


1.8 


635.529 


638.044 


0.396 


2.0 


602.282 


604.221 


0.322 


3.0 


537.444 


538.111 


0.124 


4.0 


519.496 


519.794 


0.058 


5.0 


512.080 


512.237 


0.031 



Figure 3: Graph of Xi(h) and v\(h) Table 1: Values of Ai, V\ and a 

guaranteed error estimate on \i\ 



so 



2«, = p'i {^p'i{h 2 a g )h-^j . 



Using equations (0), (0), ([21]), (0), and (0) we see that 

= cr(/i, a ff , 1) — ^dfCtg + 2a/a g + pxioif) — 2aja g 

= pi(h 2 a g )h~ 4 - p' 1 (h 2 a g )h~ 2 a g + p x (^p'i(h 2 a g )h^ . 

Inequality (|ll|) implies that 

Pi Qp'iO 2 ")^ 2 ) > 2vr 2 > 2a 



for a < n 2 and 



Pi Qp'i^ 2 ")^ 2 ) < 2^c 2 < 2a 



for a > dc 2 , so there is at least one solution of equation ([36]) in the interval 
[7T 2 ,<ic 2 ], and there are no solutions outside. Hence a g G [7r 2 ,dc 2 ]. □ 

Note 8: Numerical evidence strongly suggests that there is a unique solution 
of equation ( j36|) for any value of h. Plots of p"(a) suggest that 

-0.0603 < p'[(a) < 
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for all positive a. The fact that p"(a) is negative is clear because we showed 
that p is concave in theorem |] (i). A unique solution of equation (|36|) is 
guaranteed however under the weaker requirement that 

-2 < p'l(a) < 

for all positive a. For then 



and so 



^-(2a) = 2< 1 -p'l(^p' 1 (h 2 a)h-^ 

^(p;Qp;w)) 

p' 1 [^p' 1 {h 2 a)h-^j-2a 



(h 2 a) 



is a strictly increasing function. □ 

Figure [3] and table [I| have been created by Mathematica using formula (24) 
for Ai, and formulae (|35|), (|36|) for v\. For every value of h taken there was, 



as expected, only one solution of equation (|3q) . As in |Dav95a|| , the Hartree 
approximation / eg) g gives the worst approximation when there is an extra 
rotational symmetry, as for h — 1. 



In theorem IT3 we show that both Ai and v\ converge to c 4 ~ 500.564 and 



the guaranteed percentage error 

'v \ ~ Ai 
Ai 



x 100 



on pi is of order h 3 . This allows us to prove the asymptotic formula (0) up 
to the same order. 



5. The Groundstate 

Theorem 9: Let /-f be the negative part of the groundstate fx of A 2 |dir 
acting in L 2 (Rh) ■ Then 

ll/fl|2 <^€ (39) 



IIAlb (A 3 -^i) 1 / 2 



18 



where Xi(h), X^ih) and v\{h) are given by theorem [] and equation 

We give an asymptotic expansion of the above bound in corollary |14]. It is 
also possible to evaluate this bound for smaller values of h giving the results 
in figure || and table |[ 

Proof: There exists a complete orthonormal sequence of eigenfunctions f n 
of the biharmonic operator acting in L 2 (Rh) with corresponding eigenvalues 
\i n listed in increasing order. Let {f,g) minimize £. Then 

||/® g\\ 2 = 1 and \ x < ^ = \\/±hf 2 < ||A(/ ® g)f 2 = v x . (40) 

By writing / ® g — Y^=i a nfn with a x > 0, we see that 

oo oo 

y]al = l, a 2 = and Ai < ji\ < afy n = v x . (41) 

n=l n=l 

Using the above information, 

oo 

V\ - \l\ = ®lHn ~ (1 - «i)^l 

n=3 

oo 

> h^2<*i - (i - «i)/ii 

71=3 

= (A3-A*i)(l-a?). (42) 

Rearranging (42), we see that 

1 + «! = I — (1 — QtiJ > I — > — . 

A3 — A3 — /ii 

Using (42) again, 

1 - «i < ^^r- (43) 
2(A 3 - z/i) 

Now 

00 

ll/®^-/i||^ = (l-«i) 2 + ^^ 

n=3 

00 

= + -2^) 



n=l 
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0.05 



0.04 



. 03 



0.01 



h 



h 


(As-n) 1 / 2 


1 


0.0484 


10 


0.0336 


20 


0.0249 


40 


0.0179 


60 


0.0147 


80 


0.0128 


100 


0.0114 



Figure 4: Graph of {v x - A 1 ) 1 / 2 /(A 3 - ui) 1/2 

= 2(l-ai) 
z/i - Ai 



Table 2 



< 



A 3 - v\ 

Using theorem [l] (ii) we see that / <g> g is positive. It follows that 

l/ri<i/®s-A|. 



(44) 



□ 



The following Sobolev embedding theorem is used to convert the L 2 bound 
above to an L°° bound. 

Lemma 10: Let f G W 2 ' 2 (R 2 ). Then 



1 

< - 



i\\*fh- 



loo - 41 

Proof: Let h G L 2 (M 2 ) be the function whose Fourier transform is 

M0 = (2vr)- 1 ( 7 + 7 - 1 |e| 4 )- 1/2 . 

Then 

M = lM = g- 

Now 

( 7 + 7 - 1 A 2 )- 1 /^ = h*g 

so 



for all # G L 2 (R 2 ] 

\\(i + i- 1 ^)- 1/2 g\L = \\h*g\\l < 

20 



1 



(45) 
(46) 
(47) 
(48) 

(49) 



Hence for / e W 2 ' 2 (R N ) 

\\f\L < gll(7 + 7 -1 ^ 2 ) 1/2 /ll2 = |(7ll/ll2 + 7~ 1 ||A/|| 
The result is obtained by setting 

7 Wfh ' 

Theorem 11: Using the notation of theorem®, 



□ 



II/1II2 2(A 3 -, 1 ) 1 / 2 " (5 ° J 
Proof: Let /, g, and fx be as in theorem ^|. Then 

oo 

II A(/ ® <7 - A) 111 = a nVn + 2(1 - ai)fii 

n=l 

< z/i - Ai + ) i/j. 

VA 3 - v\ 

iyx - AO x 



(A 3 - fi) 
Lemma |1(] implies that 

||/®^-/i||^<^||/®^-/i||ll|A(/®^-A)|| 2 
lb 

iyx - Ai) 2 A 3 

< 16(A 3 -^i) 2 ' □ 



b. Asymptotic Estimates 

The functions Ai and V\ become increasingly accurate estimates of \i\ as 
h —>■ oo (see figure [3] and table [I|). The following theorem is proved using the 
fact that Ai and v\ have the same asymptotic formulae up to 0(h~ 3 ). 

Theorem 12: The first eigenvalue of the biharmonic operator acting in 
L 2 (Rh) has the asymptotic formula 

^(h) = c 4 + 2dc 2 7r 2 h- 2 + 0(h- 3 ) (51) 
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as h — > oo. 

Proof: Let a(h) be the solution of (|36| ) which makes (|35D valid. It follows 
from inequality (|TT| ) that 

^(/» a o(fc))r a = o(/r 2 ) 

as h — > oo. Hence by equation (|3T)|), 

a(/i) = rf c 2 + 0(/T 2 ) (52) 

as h — > oo. Substituting this expression into the formula (|35|) for z/i we see 
that 

Vl (h) = h^ Pl (dc 2 h 2 + 0(1)) - h- 2 p[(dc 2 h 2 + 0(l))(dc 2 + 0(h- 2 )) 

+ Pl (p' 1 (dch 2 + 0(l))/2h 2 ) 
= h-\2it 2 dc 2 h 2 + 4V2n 2 d 1/2 ch + 0(1)) 

- /T 2 (2tt 2 + 2v / 2VtT 1/ V 1 /r 1 + 0(h- 2 ))(dc 2 + 0(/T 2 )) 

+ p(7r 2 / i - 2 + v^^W 3 + 0(h- 4 )) 
= 27r 2 dc 2 h~ 2 + 4v / 2vr 2 rf 1 / 2 c/i- 3 - 2-a 2 dc 2 h- 2 - 2V2n 2 d 1 / 2 ch- 3 

+ c 4 + 2dc 2 7i 2 h- 2 + 2v / 27r 2 c/ 1 / 2 c/r 3 + 0(h- 4 ) 
= c 4 + 2dcV/T- 2 + 4 v / 2Vd 1 / 2 c/T 3 + 0(/T 4 ), (53) 

as /i — > oo. 

Also, by substituting the asymptotic formulae (10) for p\ into the formula (24) 
for Ai, we see that 

\ x (h) = c 4 + 2dc 2 ir 2 h- 2 + h- 4 (2n 2 h 2 n 2 + A^^hn + 0(1)) - 2ii A br 2 + 0(/T 4 ) 
= c 4 + 2dc 2 n 2 h- 2 + AV2n 3 h- 3 + 0(/i~ 4 ), (54) 

as h — > oo. □ 

Note 13: For long thin rectangles a good approximation to the groundstate 
of the biharmonic operator is 

V2h-V 2 sm(^) gi (y) 
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where g± is defined by formula (P3|). Intuitively one expects this function 
to be a fairly good approximation because the boundary conditions at the 
ends of the rectangle become less influential on the eigenf unction. Note that 
the above function does not actually lie in the quadratic form domain. The 
energy of this function may be computed however if we ignore this fact, and 
we see that 

Q(f) = C A + 2dc 2 7l 2 h~ 2 + 7r A h~ A . 



This compares well with the asymptotic formula (51]) for /ii as h — > oo, 
differing only by terms of order h~ 3 . □ 



Corollary 14: The bounds (39) and fl50"|) have asymptotic formulae 
IIT II' < <£ZWL = 2^(^c-,)V^_ 1/2 + 0(ft _ 3/2) 



IIAIh (As-^i) 1 / 2 27T 



ll/rll oo , (^-Ai) 1/2 Af _ 2^_^ 
IIAIh 2(A 3 -^ " 4vr 



(56) 



as h — > 00. 

Proof: The asymptotic formula of A3, 

A 3 (/i) = c 4 + (2dcV + 167r 4 )/T 2 + 0(/T 3 ) (57) 

as /i — > 00 is found by using formula (25) and the asymptotic formulae (10) for 
Pi and P3. The corollary follows by using the formulae (53), (54) and (|57|). □ 

Note 15: It is conjectured that 

IIAIh = cosh(c/2) sinh(c/2) - cos(c/2) sin(c/2) fel/2 ^-i^x 
IIAIU 2 v / 2[cosh 2 (c/2)sin(c/2) + sinh(c/2) cos 2 (c/2)] 

as /i — > 00 so 

IIA ll oo _ || A W oo IIAIh 

IIA Hoc ~ HA Ih II/: 



1 \\oo 



2 1 / 4 (d 1 / 2 c - 7r) 1 / 2 C rcosh(c/2) sinh(c/2) - cos(c/2) sin(c/2)l , 

< 1= — p 7TT ; - + 0{h 

8 v / 27r[cosh 2 (c/2) sin(c/2) + sinh(c/2) cos 2 (c/2) 1 
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< 0.121 



(59) 



for h large enough. Comparing this expression with (55) we see that bound (|5p 
is of some use, but is a lot weaker than Q59|). An improved bound would be 
desirable. □ 



7. Appendix 



Proof of theorem [j] (v): 

Preliminary calculations show that the roots (3 n of equation (Q) are of the 
form a + n 2 TT 2 + o(l) as a — > oo. Define 

P + (a) = a + n 2 n 2 + 2 v / 2n 2 7r 2 a- 1 / 2 + GnVa' 1 + -(-l) n 5V2n 4 n 2 a- 3 / 2 . 

6 



Then 



cos 



a 



a: 



VP. 



2 - a 2 



tanh \/ P + + a sin \[~fi^ 



a: 



cos(n 2 vr 2 + 2v / 2Vvr 2 C r 1 / 2 + O^- 1 )) 1 ' 2 



1 + 2v / 2V 1/2 



12 + n 2 7r 



2wk 

x tanh(2a + 0(l)) 1/2 



-a 



+ 0(a" 3 ^ 



-1/2 



x sin |n 2 7r 2 + 2v / 2Vvr 2 a~ 1/2 + §n 2 -n 2 cT x + ^(-l) n 5v / 2' 
(-^'(l-AV + Ofa- 3 / 2 )) 

_ V2^! (1 _ ^-1/2 _ l 2f 2 a -i + (a- 3 / 2 ))(l + 0(a- 3 / 2 )) 



n 4 7r 4 a- 3 / 2 



1/2 



2ri7r 



4 



x <-l) n ( V2mra- 1/2 + 2wroT 1 - 2v / 2mra~ 3/2 - -v^nVcT 3 / 2 



^(-^"Sv^nVa-^ + O^- 2 ) 



-ir 



12 

< sech v/3+ + a 



4 5nV(l 



a 1 + 0(« 



-3/2i 
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for a large enough. 
Similarly, defining 



fi_(a) = a + nV + 2y/2: 



'n 2 Tr 2 a~ 1/2 + Qn 2 7t 2 a 



-(-l) n 5V2n 4 7r 4 a 



-3/2 



we see that 



cos sffi. 



a 



tanh + ct sin — at 



L. - a - 



(-l) n 4 



_ J_ 5n V(l-(-l) n ) a^ + 0(a- 3 / 2 ) 



> sech + a 



for a large enough. 

It follows that for a large enough, (3 n lies between /3_ and (3+ so 



/9 n (a) = a + nV + 2v / 2nVa~ 1/2 + Qn 2 n 2 a~ l + 0(a~ 3/2 ). (60) 

Thus 



Let F(ot,P) be the left hand side of equation @. Then differentiating @, 
we see that 



2(F 1 (a, AQ + F 2 (q, &))(/% - a 2 ) 3 / 2 
cosh y/ p n + a COS yj (3 n — a 
= 2(p n - a) {f3 n + a) 1/2 tanh 

— 2a(f3 n — a) 1//2 tan \J (5 n — a 

- 2(3 n ((3 n - a) l/2 ((3 n + a)~ 1/2 tanh v^ + cttan v 7 ^ - a 



p n (a) = (3 n (a) 2 - a 2 

= 2n 2 n 2 a + 4v / 2n 2 vrV /2 + 0(1). 
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... = _2n 2 7r 2 + 0(a- 1 / 2 ). 

Also, 

2F ^f nWn ~77p- = - 2 «)(^ - «) 1/2 (^ + «) tanh V^+^ 

COSh A/Pn + « COS VPn — a 

- (p n + 2a)(P n - a)((3 n + a) 1/2 tan y/J n - a 
+ 2a(3 n tanh \ j (3 n + a tan (3 n — a 
... = ->/2a 3/2 + 0(a). 

Hence using (pip 

^ = 1 " I -^/ 2 + (a) 

= l-v / 2n 2 vr 2 ^ 3 / 2 + 0(a- 2 ). (62) 



Therefore 



p' n (a) = 2(3 n (3' n - 2a 

= 2n 2 n 2 + 2V2n 2 n 2 a~ l/2 + 0(a~ 1 ). □ 
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